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Interacting fermions are ubiquitous in nature and understanding their ther- 
modynamics is an important problem. We measure the equation of state of a 
two-component ultracold Fermi gas for a wide range of interaction strengths at 
low temperature. A detailed comparison with theories including Monte-Carlo 
calculations and the Lee-Huang- Yang corrections for low-density bosonic and 
fermionic superfluids is presented. The low-temperature phase diagram of the 
spin imbalanced gas reveals Fermi liquid behavior of the partially polarized 
normal phase for all but the weakest interactions. Our results provide a bench- 
mark for many-body theories and are relevant to other fermionic systems such 
as the crust of neutron stars. 
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Recently, ultracold atomic Fermi gases have become a tool of choice to study strongly cor- 
related quantum systems because of their high controllability, purity and tunability of interac- 
tions (7). In the zero-range limit, interactions in a degenerate Fermi system with two spin- 
components are completely characterized by a single parameter l/kpa, where a is the s-wave 
scattering length and kp = (Gvr^n)^/^ is the Fermi momentum (n is the density per spin state). 
In cold atom gases the value of |a| can be tuned over several orders of magnitude using a Fesh- 
bach resonance, this offers an opportunity to entirely explore the so-called BCS-BEC crossover, 
i.e. the smooth transition from Bardeen-Cooper-Schrieffer (BCS) superfluidity at small negative 
values of a to molecular Bose-Einstein Condensation (BEC) at small positive values of a {1, 2). 
Between these two well-understood limiting situations a diverges, leading to strong quantum 
correlations. The description of this system is a challenge for many -body theories, as testified by 
the large amount of work in recent years (i). The physics of the BEC-BCS crossover is relevant 
for very different systems, ranging from neutron stars to heavy nuclei and superconductors. 

In the grand-canonical ensemble and at zero temperature, dimensional analysis shows that 
the Equation of State (EoS) of a two-component Fermi gas, relating the pressure P to the chem- 
ical potentials jii and 112 of the spin components can be written as 



where Po(a*i) = l/lf>n'^{'2m/h^fl'^ixT is the pressure of a single-component ideal Fermi gas, 
m is the atom mass and fi is the Planck constant divided by 27r. The indices 1 and 2 refer to the 
majority and minority spin components, respectively. From the dimensionless function h{5i, rf), 
it is possible to deduce all the thermodynamic properties of the gas, such as the compressibility, 
magnetization or the existence of phase transitions; the aim of this paper is to measure 77) 
for a range of interactions and spin imbalances {rf) and discuss its physical content. 5i is 
the grand-canonical analog of the dimensionless interaction parameter l/kpa. 




(1) 



2 



In-situ absorption images of harmonically trapped gases are particularly suited to investigate 
their EoS as first demonstrated at MIT {3) and ENS (4). In the particular case of the grand- 
canonical ensemble, a simple formula relates the local pressure P at a distance z from the 
center of the trap along the z axis to the doubly-integrated density profiles ni and n2 (5): 
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P{f,,{z),f,2{z),a) = {Mz)+M^)) ■ (2) 

Here we define the local chemical potentials ijLi{z) — /i° — \mulz'^, where ji^ is the chemical 
potential of the component i at the bottom of the trap, assuming local density approximation. 
Uj. and are the transverse and axial angular frequencies of a cylindrically symmetric trap 
respectively, and ni{z) — J ni{x,y, z)6xdy, is the atomic density rii of the component i, 
doubly integrated over the transverse x and y directions. In a single experimental run at a given 
magnetic field, two images are recorded, providing ni{z) and n2{z) (see Fig.S4 in (6)); the 
z-dependence of the chemical potentials then enables the measurement of P along a curve in 
the Tj) plane (6). This method was validated in (4) for the particular case of the unitary limit 
a — oo. Deducing the function h from the doubly integrated profiles further requires a precise 
calibration of and the knowledge of the central chemical potentials /i° (6). 

Our experimental setup is presented in (7). We prepared an imbalanced mixture of ^Li in 
the two lowest internal spin states, at the magnetic field of 834 G (where a — oo), and trapped 
it in a hybrid magnetic-optical dipole trap. We then performed evaporative cooling by lowering 
the optical trap power, while the magnetic field was ramped to the final desired value for a. 
The cloud typically contained A/^ = 2tol0xl0'^ atoms in each spin state at a temperature of 
0.03(3) Tp, justifying our T — assumption (6). The final trap frequencies are ujz/27r ~ 30 Hz, 
u!r/27r ~ 1 kHz. Below a critical spin population imbalance, our atomic sample consists of a 
fully-paired superfluid occupying the center of the trap, surrounded by a normal mixed phase 
and an outer rim of an ideal gas of majority component atoms (4, 7, 8). 
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For a given magnetic field, 10 to 20 images are taken, leading after averaging to a low- 
noise EoS along one line in the {Si,r]) plane. Measurements at different magnetic fields chosen 
between 766 G and 981 G give a sampling of the surface h(5i, r]) in the range —\<5i< 0.6 
and —2<r]< 0.7 (Fig{T]). Let A(5i) be the limiting value of the ratio of chemical potentials 
fii{z) / fi2{z) below which the minority density vanishes. At fixed 5i and r] < A(5i), h(5i,r]) 
represents the EoS of an ideal Fermi gas of majority atoms and is equal to 1. For i] > A{Si), 
it slowly rises and corresponds to the normal mixed phase, where both spin components are 
present. At a critical value i] = ?7c(5i), the slope of h abruptly changes (6), the signature of a 
first-order phase transition from the normal phase (for A < r] < rjc) to a superfluid phase with a 
lower chemical potential imbalance (r/ > r]c). We notice that the discontinuity is present for all 
values of 5i we investigated, and this feature is more pronounced on the BEC side. 

Let us first consider the EoS of the superfluid phase, rj > 7]^. Each of our in-situ images 
has, along the 2-axis, values of the chemical potential ratio ri{z) = H2{z) / ni{z) both lower 
and greater than r]c. In the region where r](z) > rjc the doubly-integrated density difference 
ni{z) — n2{z) is constant within our signal-to-noise ratio (see Fig.S4). This is the signature of 
equal densities of the two species in the superfluid core, i.e. the superfluid is fully paired. Using 
Gibbs-Duhem relation rii = equal densities rii = n2 imply that P(/ii, ^2, a) is a function 
of /i and a only, where /i = (/ii + /i2)/2. For the balanced superfluid, we then write the EoS 
symmetrically: 



In order to avoid using negative chemical potentials, we define here /i = jj, — Eh/2, where 



then a single- variable function. It fully describes the ground state macroscopic properties of the 
balanced superfluid in the BEC-BCS crossover and is displayed in Fig|2]as black dots. 

In order to extract relevant physical quantities, such as beyond mean-field corrections, it is 
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(3) 



is the molecular binding energy Ei, 




2 for a > (and for a < 0). hs{5) is 



convenient to parametrize our data with analytic functions. In this pursuit, we use Pade-type 
approximants (6), interpolating between the EoS measured around unitarity and the well-known 
mean-field expansions on the BEC and BCS limits. The two analytic functions, h^'^^ and h^^'~^ 
are respectively represented in blue and red solid lines in Fig l2] and represent our best estimate 
of the EoS in the whole BEC-BCS crossover. 

First, on the BCS side 5 < 0, h^*^^ yields the following perturbative expansion of the energy 
in series of kpa: 



where is the total number of atoms, Ep is the Fermi energy and where by construction of 
h^'^^, the mean-field term (proportional to kpo) is fixed to its exact value IO/Qtt. We obtain 
beyond mean-field corrections up to 3^'^ order. The term proportional to {kpaY agrees with 
the Lee- Yang (9, 10) theoretical calculation 4(11 - 21og2)/2l7r2 ~ 0.186. The third order 
coefficient also agrees with the value 0.030 computed in {11). 
Second, around unitarity the EoS expands as 



We find the universal parameter of the unitary T = superfluid, = 0.41(1) with 2 % 
accuracy. This value is in agreement with recent calculations and measurements (1). Our 
thermodynamic measurement = 0.93(5) can be compared with a recent experimental value 
C = 0.91(4) (12) as well as the theoretical value ( = 0.95 (13), both of them obtained through 
the study of the pair correlation function. This experimental agreement confirms the remarkable 
link between the macroscopic thermodynamic properties and the microscopic short-range pair 
correlations, as shown theoretically in (14). 

Third, in the BEC limit the EoS of the superfluid is that of a weakly interacting Bose- 





(4) 
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Einstein condensate of molecules (9, 15): 

where add = 0.6a is the dimer-dimer scattering length (i) and n is the dimer density. The term 
in \Jna\^ is the well-known Lee-Huang- Yang (LHY) correction to the mean-field interaction 
between molecules (9, 15). Signatures of beyond mean-field effects were previously observed 
through a pioneering study of collective modes {16) and density profile analysis (77) but no 
quantitative comparison with ([5]) was made. Fitting our data in the deep BEC regime with 
Eq.®, we measure the bosonic LHY coefficient 4.4(5), in agreement with the exact value 
128/15^71" — 4.81 calculated for elementary bosons in (9) and recently for composite bosons 
in (75). 

Having checked this important beyond mean-field contribution, we can go one step fur- 
ther in the expansion. The analogy with point-like bosons suggests to write the next term as 

{6, 18, 19). Using h^^^{6) (Fig|2l and (6)), we deduce the 
effective three-body parameter for composite bosons B = 7(1). Interestingly, this value is close 
to the bosonic hard-sphere calculation B = 8.5 (20) and to the value B ^ 7.2 for point-like 
bosons with large scattering length (79). 

Our measurements also allow direct comparison with advanced many-body theories devel- 
oped for homogeneous gases in the strongly correlated regime. As displayed in FigjSK, our 
data are in agreement with a Nozieres-Schmitt-Rink approximation (27) but shows significant 
differences from a quantum Monte-Carlo calculation (22) and a diagrammatic approach (25). 
The measured EoS strongly disfavors the prediction of BCS mean-field theory. 

Comparison with Fixed-Node Monte-Carlo theories requires the calculation of the EoS 
^(l/kpa) in the canonical ensemble: 



l{A7r-3V3)nald{log{nald) + B) 



that is deduced from hg^^iS) and hg^'-^iS) (6). As shown in Figl3^, the agreement with theories 
(24-26) is very good. 

We now discuss the EoS of the partially polarized normal phase (black points in Fig{T]). At 
low concentrations, we expect the minority atoms to behave as non-interacting quasi-particles, 
the fermionic polarons (27). The polarons are dressed by the majority Fermi sea through a 
renormalized chemical potential /i2 — A(5i)iii (28) and an effective mass m*{Si) (26, 29, 30). 
Following a Fermi liquid picture, we propose to express the gas pressure as the sum of the Fermi 
pressure of the bare majority atoms and of the polarons (4): 

= 1 + l^^Y' iv - Am'/'. (6) 
\ m J 

Our measured EoS agrees with this model at unitarity and on the BEC side of the resonance 
(Fig{T]), where we use for m*{5i) the most advanced calculations (30, 31). On the BCS side of 
the resonance however, we observe at large minority concentrations an intriguing deviation to 

In the BCS regime, the superfluid is less robust to spin imbalance. Consequently, the ratio of 
the two densities n 1/^2 in the normal phase becomes close to unity near the superfluid/normal 
boundary r]c. The polaron ideal gas picture then fails. 

Alternatively, we can let the effective mass m* be a free parameter in model © in the fit 
of our data around t] = A. We obtain the value of the polaron effective mass in the BEC-BCS 
crossover (FiglD). 

An important consistency check of our study is provided by the comparison between our 
direct measurements of i]c{Si) (from Fig{TJ black dots in the inset of FigjH) and a calculated 
r]c{Si) from Eq.© and the EoS of the superfluid hs. Assuming negligible surface tension, 
the normal/superfluid boundary is given by equating the pressure and chemical potential in 
the two phases. This procedure leads to the solid red line in the inset of FiglH in excellent 
agreement with the direct measurements. In addition, by integrating our measured EoS of the 
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homogeneous gas over the trap, one retrieves the critical polarization for superfluidity of a 
trapped gas, in agreement with most previous measurements (6). 

We have measured the equation of state of a two-component Fermi gas at zero temperature 
in the BEC-BCS crossover. As a first extension, we could explore the thermodynamics of the 
far BEC region of the phase diagram where a new phase associated with a polarized superfluid 
appears (17, 26). Another interesting direction would be the mapping of the EoS as a function 
of temperature and investigate the influence of finite effective range that is playing a key role in 
higher density parts of neutron stars. 
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Figure 1: h{Si, r]) of a zero-temperature two-component Fermi gas in the BEC-BCS crossover. 
(A): Samples of the data for different magnetic fields. The black (red) data points correspond to 
the normal (superfluid) phase, and are separated at r]c{5i) by a clear kink in the local slope of h. 
Solid black lines are the predictions of the polaron ideal gas model, Eq Q. The scattering length 
corresponding to each curve is (from left to right) : (1.7, 3.4, oo, —1.3) in units of 10"^ ao, where 
ttQ is the Bohr radius. (B): h{5i, r]). The black dots are data recorded for each magnetic field 
value (as in Fig{T^). The black lines correspond to the parametric curves (Si{r]),r]) scanned by 
the density inhomogeneity in the harmonic trap (6). The red line is A{5i), the frontier between 
the fully polarized (FP) ideal gas h = 1 and the normal partially polarized (PP) phase. The 
green line is 'qd^i) marking the phase transition between the normal and superfluid (S) phases. 
The surface is the parametrization of h{6i, rj) given in the text. 
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Figure 2: hs{5) of the T = balanced superfluid in the BEC-BCS crossover (black dots). The 
blue solid line is the fit hg'"^{5) on the BCS side of the resonance, the red solid line is the 
fit hg^^{5) on the EEC side (see text). The dotted (dashed) red line is the mean-field (Lee- 
Huang- Yang) theory (33). Inset: Zoom on the BCS side. The dotted (resp. dashed) blue line 
is the EoS including the mean-field (resp. Lee- Yang) term. The systematic uncertainties on the 
X and y-axis are about 5 %. The errors bars represent the standard deviation of the statistical 
uncertainty. 
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(A) (B) 




Figure 3: Comparison with many-body theories. (A): Direct comparison of hs{5) with a 
quantum Monte-Carlo calculation, red open circles (22), a diagrammatic method, green open 
squares (23), a Nozieres-Schmitt-Rink approximation, blue open triangles (21) and the BCS 
mean-field theory, solid blue line. Inset: Zoom on the BCS side. (B) EoS in the canonical en- 
semble ^(1/ kpo) (solid black line) deduced from the Pade-type approximants to the experimen- 
tal data h^'^^ and plotted in FiglH Fixed-Node Monte-Carlo theories: red squares (24), 
blue circles (25), green triangles (26). 
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-1.0 -0.5 0.0 0.5 

Figure 4: Effective mass m* /m of the polaron in the BEC-BCS crossover (black dots). The 
blue dashed line is a calculation from (29), red open squares {30), green dot-dashed line (2(5), 
and blue solid line {31). Measurements at unitarity through density profile analysis (blue trian- 
gle {3)) and collective modes study (brown empty circle (7)) are also displayed. Inset: Phase 
diagram of a zero-temperature imbalanced Fermi gas in the BEC-BCS crossover. The blue line 
is the theoretical value of A {26, 29, 30) that sets the separation between the partially polarized 
(PP) and the fully polarized (FP) phases. Black dots are the measured values of r]c (as in Fig{TK) 
which set the separation between the superfluid (S) phase and the partially polarized phase. The 
red line is the calculation of rjc using our EoS of the superfluid and the model @ for the nor- 
mal phase. The green squares are lower bounds of r/c given by the values of the gap measured 
in {32), see {6). 
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Supporting Online Material for 

The Equation of State of a Low-Temperature Fermi Gas with l\inable In- 
teractions 

N. Navon, S. Nascimbene, F. Chevy, C. Salomon 
Materials and Methods 

Extracting the EoS of the uniform gas from the density profiles of trapped clouds 

In this section, we provide additional insight on the reconstruction of the EoS of the uniform 
gas. Several steps are required to deduce the EoS from the doubly-integrated profiles. First, the 
determination of 77) along the z axis requires a precise calibration of a;^ and the knowledge 
of the central chemical potentials The axial confinement is produced by a magnetic field 
curvature, which ensures very good reproducibility, ujz is calibrated (to < 1%) by measuring the 
frequency of the axial center of mass dipole mode. \x\ is determined using the fully polarized 
outer rim of the cloud. In this region the density profile is fitted by a Thomas-Fermi formula 
n\{z) — a{\ — /R\Y'/'^, which gives ii\ — \mujlR\. h is then directly obtained by taking the 
ratio (ni(z) + n2{z))/a{l — z^ /R\Y/'^, thus avoiding the measurement of the radial frequency 
Ur and cancelling many systematic effects such as imperfect atom counting (i, 4). 

The determination of requires some information on the EoS. We use the outer radius of 
the minority component as a reference. Indeed, in the limit of vanishing minority spin density, 
the minority chemical potential is equal to the one of a single minority atom immersed in a 
Fermi sea of majority atoms, the so-called polaron problem. All advanced calculations (26, 29, 
30) and the measurement in (28) agree on the value of the chemical potential ratio A{5i) = 
^2/ A*i when n2 0, which is plotted in the inset of Fig.4. fi% is then fitted on each image so 
that T] — iJi2zl IJ-iz is equal to A{5i {niz)) at the point z where the minority density vanishes. 
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Within the local density approximation (LDA), the local chemical potentials along the z- 

axis vary as \iiz = jj,^ — ^muolz^ (for species i). The local interaction parameter, defined as 
^iz — fi/ y/2miiiza (where 1 is the majority component), also varies along the cloud, as well as 
the local chemical potential imbalance 77^ = ijl2z/ A*iz- Substituting z in 5i in favor of 77, we find: 



the center of the cloud and we dropped the z subscript for clarity. Each trapped density profile 
thus gives the EoS along the parametric curve {5i{r]),r}). 

The images used to reconstruct the EoS at a given magnetic field do not perfectly belong 
to the same curve {61(1]), r]). Here we quantify the systematic error produced by overlapping 
(and then averaging) the various images. To do so, we have simulated density profiles using 
our measured EoS with the distribution of initial parameters rj'j)} (where j is the image 
index) of the images used and have reproduced the reconstruction process. We then compare 
this result to the expected EoS corresponding to the mean values of ^°)}. The difference 
is less than 3 %. Moreover, the determination of using A{6i) as a reference (as explained 
in the text), leads to an additional systematic error. We estimate it to be 4 % on the hs axis of 
Fig.2. The uncertainty on the imaging system magnification leads to a 5 % systematic error on 
the 6 axis of Fig.2, while trap anharmonicity effects are expected to be 2%. 

The critical chemical potential rjc is extracted from each equation of state h{5i,r]) obtained 
at a given magnetic field. We fit the data in a region rj E [fj — 0.2, fj + 0.2] with a continuous 
function made of two straight segments, and observe that the location of the breaking point is 
insensitive to the fj value defining the set of points chosen for the fit. This shows that our data 
supports an abrupt change of slope, and rjc is identified as the breaking point. 




(7) 



where 5° = h/J2miXia (resp. rjo — 11%/ //?) is the interaction strength (resp. local imbalance) at 
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Parametrization of the Superfluid Equation of State 

In order to extract physical parameters in the BEC-BCS crossover from our equation of state 
and to calculate the canonical EoS, we use a simple parametrization of our data that possesses 
the correct asymptotic behaviors (at 5 — )■ ±00 and 5 — )■ 0). 

First, on the BCS side 5 < 0, we use a Pade-type approximant: 

hTHs) = ^~ (8) 

52 + as6 + a4 

Using the mean-field asymptotic behavior hs{6) ~ 1 — 5/ {3n5) in the BCS regime as a con- 
straint on the ai, a fit of our data for 5 < 0.2 with ([8]) leads to the coefficients gathered in 
Table S3. 

Second, on the BEC side 5 > 0, we capture the behavior in the BEC limit (Eq.(5) in the 
paper with the additional log term) using the following formula: 

^BEc.^^ ^ /9i + f32S + 1336 logjl + 6) + f3,P + f3,6^ 

The /3i being contrained by the values of and ( previously determined from the BCS side 
and by the exactly known coefficients in (5), we fit our data for 5 > —0.2 with a single free 
parameter in ^ and obtain the values of /3i in Table S3. 



The three-body B parameter 

For a weakly interacting Bose-Einstein condensate, the ground state energy can be written as 
(18,19): 



E = N- 



m 



-n\l + 



128 

15v^ 



yna? + 



-(47r - 3V3)na^(log(raa^) + B) 
3 



(10) 



The first term is the mean-field contribution, the second is the Lee-Huang- Yang correction. The 
third term involves the three-body problem and was first calculated in {18). It was shown in {15) 
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that this equation of state is rigorously vaUd for composite bosons up to the LHY term by sub- 
stracting to the energy the binding energy of the molecules, replacing a by the dimer-dimer 
scattering length qm, tti by the dimer mass and considering n as the dimer density. The anal- 
ogy with point-like bosons suggests to write the next term for a Bose-Einstein condensate of 



the low-momenta physics (18) where the composite nature of the dimers should not play a role. 
The coefficient B also involves the three-boson problem at high momenta, unveiling the inner 
structure of the dimers. For elementary bosons, B depends on the microscopic details of the 
interaction potential between bosons and involves Efimov physics (19). On the other hand, the 
internal structure of the dimers is completely characterized by the scattering length a. Hence, 
we could expect the three-dimer problem to be solely described by a as well. B would take a 
universal value, i.e. independent on the fermionic species. Using hg^'~^{d), we deduce the ef- 
fective three-body parameter for composite bosons B — 7(1). Interestingly, this value is close 
to the bosonic hard-sphere calculation B — S.b (20) and to the value B 7.2 for point-like 
bosons with large scattering length (19). 

Derivation of the pressure formula (2) 

We consider a mixture of species i, of mass m^, trapped in a harmonic trap of transverse fre- 
quencies Uri. Using Gibbs-Duhem relation at a constant temperature T, dP — J2i f^idjii, then 



where we have used local density approximation (//^(r) — /jl^ — V{r)) to convert the integral 
over space to an integral on the chemical potentials. The integral is straightforward and yields 



dimers as |(47r — 3\/3)na|^(log(na|^) + B) . The coefficient before the log term is given by 




(11) 
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Eq (2) of the main text is a special case of Eq (fTTI) for two spin components of equal masses 
confined by the same trapping potential. 



Critical Polarization for a Trapped Gas 

Here we compare our measurements of the chemical potential imbalance for the normal/superfluid 
transition to previous works (SI, S2, 7, 8), where this phase transition was characterized by mea- 
suring the maximum polarization Pc above which the superfluid is no longer present. We define 
the polarization P = {Ni — N2)/ {Ni + N2), where Ni is the total atom number for the species 
i. Modeling the normal phase using 



we calculate the atom number Ni by integrating the densities rii = dP/dni over the trap: 



The critical polarization is obtained when the transition normal to superfluid occurs at the bot- 



tom of the trap, i.e. for /i^//^? = ^ci^i), where 5° = h/ ^2m^ia. rjc{5i) is provided by 

the solid red line in Fig. 3 (see text). In Fig. SI, we plot P^. as a function of 1/kpa, where 
kp = h{ujzCo'^Y^^{6NiY/^, providing a direct comparison with previous experimental data 
(SI, S2, 7, 8). Our results are in excellent agreement with (SI, 7, 8) but not with (5*2), where 
the partially polarized phase is absent. As the atom numbers and trap anisotropy in our experi- 
ment are close to those in (52), this discrepancy remains to be understood. 

Grand Canonical-Canonical Correspondence 

In this section, we explicit the conversion formulas between the grand-canonical and canonical 
ensembles. We recall that the energy density S = E/V is linked to the pressure by the usual 




(12) 
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Legendre transform 8 = —P + nn where n = a^^tm the chemical potential in the fully paired 
superfluid and n is the total atomic density. Combining this relationship with the Gibbs-Duhem 
relation dP = ndfi and using the proper normalization for the dimensionless function h{5), one 
finds the two formulas: 

h{6)-lh'{6) 

- ^h{6) - lh'{6)r/^ ^^^^ 
Eq (fT3]) relates the natural variables of both ensembles, while Eq (fT4l) provides the canonical 
EoS. The canonical EoS (displayed in Fig. 3b) is then simply obtained through a parametric plot 
of {x{5),^{5)) using the Pade-type fitting functions for h{S). Another familiar quantity that can 
be calculated from our fits is the chemical potential in the BEC-BCS crossover: 

^ = (15) 

where 6 (6) is the Heaviside step function. 



Thermometry 

We perform thermometry of our imbalanced gases on the fully polarized outer shell of the cloud, 
which is non-interacting. We thus fit the wings of the density profiles with finite-temperature 
Thomas-Fermi distributions (53). Using this technique, we find a temperature of T/Tp = 
0.03(3) at the unitary limit and on the BEC side of the resonance, a temperature much smaller 
than the critical temperature for superfluidity in this interaction range, justifying the T = 
assumption. 

On the BCS side, the small size of the fully polarized region renders this method inaccurate 
and gives an experimental upper bound T/Tp < 0.13. However, the observation of a wide 
superfluid plateau on the doubly-integrated density difference is a signature of an unpolarized 
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superfluid phase clearly indicating that T < Tc.ln the normal phase, finite temperature correc- 
tions are expected to be of the order of {T/Tpf. As T < Tc < Tp in the BCS regime, they 
are expected to be negligible. In addition, we observe an abrupt change of slope at the nor- 
mal/superfluid boundary indicating the presence of a first-order phase transition. Consequently, 
this gives an upper bound T < Ttri where Ttii is the temperature of the tri-critical point at which 
the phase transition becomes of second order. Finite-temperature corrections are expected to 
be dominated by thermal excitations of Bogoliubov-Anderson phonons. From the superfluid 
equation of state determined from our data we compute the speed of sound c = ^n/mdii/dn 
as a function of interaction strength, and the finite-temperature correction to the pressure for 
T — Ttri predicted by a mean-field theory {S4, S5, S6). We thus infer that the systematic error 
on hs due to a finite temperature is less than 2 % for our data on the BCS side of the resonance. 

Link between the critical imbalance rjc and the pairing gap A 

An intriguing link exists between the measurements of ric{6i) and the superfluid pairing gap 
A. Indeed, for Si < 0.6 we have observed that the superfluid is unpolarized. An unpolarized 
superfluid becomes unstable as soon as flipping the spin of a minority atom decreases the grand- 
potential Q. = E — HiNi — 112^2- After the spin-flip, a pair is broken and releases an energy 
equals to A£' = 2A. Thus the total grand-potential difference is AQ = //2 — A*i + 2A. The 
superfluid is stable against infinitesimal polarization as long as AQ > 0, hence the minority 
chemical potential jjL2 must be lower than fi — A, where fi = {{xi + )U2)/2. This argument yields 
lower bounds on the values of r^c^ Vc > (1~^/a*)/(1 + ^/a*) (empty green squares in the inset 
of Fig.4). We observe that our rjc data is very close to these bounds, calculated from previous 
experimental determinations of the gap from (32). As initially pointed out in (S7) for the unitary 
gas, this suggests that in the investigated range the superfluid to normal transition is driven by 
single-particle excitations of the superfluid. 
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\lkpa 

Fig. SI. Critical Polarization Pc of a trapped imbalanced Fermi gas in the BEC-BCS crossover 
(solid black line) and comparison with the data from MIT (red squares (S) and blue triangles 
{SI)), ENS (green circle (7)), and Rice (empty black square (52)). In the gray region P > Pc, 
the superfluid phase is absent. 
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a2 


as 


0:4 






-1.137 


0.533 


-0.606 


0.141 
















3.78 


8.22 


8.22 


-4.21 


3.65 


0.186 



Table S2. Pade-type approximants coefficients ctj and Pi fitted from our data. 
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Fig. S3. Canonical-Grand Canonical Correspondence in the BEC-BCS crossover, (a) Canonical 
natural variable l/kpa as a function of the grand-canonical natural variable 5. (b) Chemical 
potential in the canonical ensemble (black solid line). For comparison, mean-field BCS theory 
is the dashed blue line, the binding energy, the dashed red line. 
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z 

Fig. S4. Raw data: doubly-integrated density profiles n2{z) (gray dots) and nd{z) = ni{z) — 
n2iz) (black dots) for a gas prepared in the unitary limit. Ri, R2, Rs are the boundaries of the 
fully polarized phase, of the partially polarized phase, and of the superfluid core, respectively. 
The plateau on the density difference nd{z) observed in the region \z\ < Rs indicates equal 
densities for both spin components in the superfluid phase. 
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